RD-A124  970  IDENTIFICATION  D' UNE  CLASSE  DE  PROCESSUS  DE  POISSON  ' 
FILTRES  (IDENTIFICRTI.  .  (U)  NORTH  CAROLINA  UNIV  AT 
CHAPEL  HILL  DEPT  OF  STATISTICS  D  D  0RUCQ  20  MOV  82 
UNCLASSIFIED  N00014-75-C-0491  F/G  1271 


Ab.AtfV  97 o 


IDENTIFICATION  OF  A  CLASS  OF  FILTERED  POISSON  PROCESSES 
(Identification  d'une  Classe  de  processus  de  Poisson  filtres) 


1 
l 


C; 


Denis  de  Brucq** 


‘Research  supported  by  ONR  Contracts  N00014-75-C-0491  and 
N00014-81-K-0373. 


* ‘University  of  Rouen,  Laboratory  for  Signal  Processing,  BP  67,  76130 
Mt.  St.  Aignan,  France  and  University  of  North  Carolina,  Department 
of  Statistics,  Chapel  Hill,  NC  27514. 


±4 

« 

1 

1 

1 

#  •< 

fii 

c-: 


CL. 

O 

*  .  > 

C_5 

St'  : 

UJ 

X  •  A8  1983  ’  M 

_ 1 

U. 

&  Jr 

g 

A 

dulnbution  is  unlimi.cvd 

83  11  28  055 

] 


■/( 


SECURITY  CLASSIFICATION  of  this  PAGE  /Phan  Dote  Untered) 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

1.  RECIPIENT'S  CATALOG  NUMBER 

4.  TITLE  (end  Subtitle) 

^Identification  of  a  Class  of  Filtered  Poisson 
Processes  (Identification  d'une  Classe  de 
processus  de  Poisson  filtres) 

S.  TYPE  OF  REPORT  «  PERIOD  COVEREO 

Technical 

S.  PERFORMING  ORG.  REPORT  NUMBER 

7.  author^; 

Denis  de^Brucq 

a.  contract  or  grant  NUMBERr*; 

N00014-75-C-0491 

N00014-81-K-0373 

».  PERFORMING  ORGANIZATION  NAME  AND  AOORESS 

Department  of  Statistics 

University  of  North  Carolina 

Chat>el  Hill.  North  CaroTina  77514 

10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  S  WORK  UNIT  NUMBERS 

NR  042-269 

SRO  105 

n.  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Statistics  6  Probability  Program 

Office  of  Naval  Research 

Arlington.  VA  22217  - 

12.  REPORT  DATE 

13  NUMBER  OF  PAGES 

7 

14.  MOMITQRIMG  AGENCY  name  A  AOORESS/i/  dl Iterant  from  Controlling  Olltce) 

15.  SECURITY  CLASS,  (ol  thle  report) 

UNCLASSIFIED 

ISa.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 

t«.  DISTRIBUTION  STATEMENT  (ol  thle  Report) 

Approved  for  public  release;  distribution  unlimited 

17.  DISTRIBUTION  STATEMENT  (ol  tlte  ebetract  entered  In  Ulock  JO.  II  dl Iterant  Irum  Report) 

u.  supplementary  notes 

II.  KEY  WORDS  (Continue  on  r««r«r««  it  n«e«a««rr  ond  identity  by  block  number) 

Filtered  Poisson  processes;  data  modeling  and  identification,  underwater 
acoustics. 

20.  ABSTRACT  (Continue  an  reeeree  elde  II  neceeeary  and  Identity  by  block  number) 


^  The  formulation  of  a  model  is  an  essential  step  in  the  experimental 
process  of  understanding  a  phenomenon.  In  this  paper  a  new  class  of 
filtered  Poisson  processes  is  introduced:  the  amplitude  has  a  low  which  is 
spherically  invariant  and  the  filter  is  real,  linear  and  causal.  It  is  shown 
how  such  a  model  can  be  identified  from  experimental  data. 


IDENTIFICATION  o'UhE  cuasse  de  PROCESSUS  DE  POISSON 
FILTRES 


Deni*  oe  BRUCQ 


University  oe  ROUEN,  Laboratoire  oe  Traitemenl  oe  ('Information,  BP67  76130  MONT  SAINT  AIGNAN 


RESUME 


SUMMARY 


La  formulation  o'un  moo4le  eat  un  pd*  easentiel  pour  The  formulation  of  a  mooel  I*  an  essential  step  in 
la  comprehension  o'un  phenomena  experimental.  Dans  the  experimental  process  of  unoerstanolng  a  pheno- 
ce  travail,  une  nouveii*  class*  oe  processus  oe  Pols-  menon.  In  this  paper  a  new  class  of  filtrea  Poisson 
son,  filtr4s  est  introtxilt*  :  I'amplituae  aXatoire  suit  processes  Is  Introoucea  :  the  amplltuoe  has  a  law 
une  loi  spheriquement  invariant*  *t  la  flltre  est  reel,  which  Is  spherically  Invariant  ana  the  filter  is  reel 

lineaire,  causal  4  minimum  am  phase.  Avec  ces  hypo-  linear  ana  causal.  It  Is  shown  how  such  a  model  can 

theses,  le  mooeie  peut  etre  loentifie  4  partlr  oes  r6-  be  laentifieo  from  experimental  oata. 
sultats  exp4r imentaux. 
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Introouctlon 

Le  module  introouit  cans  ce  travail  est  un  outil  uti¬ 
le  pour  o6crire  par  example  le  bruit  acoustique  marln, 
le  bruit  o'appareillages  eiectroniques,  certalnes  par¬ 
ties  o'tflectroenclphalogrammes. 

Un  phenomena  initial,  consloere  comma  poissonnien 
oans  le  temps,  a  une  amplituoe  altfatolre  C  oe  loi 
sph4r iquement  Invariants.  Ce  phenomena  est  observe 
au  moyen  o'un  filtre  linlaire  causal  oe  r<ponse  Impul- 

sionnelle  G.  Ainsi  ('observation  preno  la  forme 

+• 

7<ui,  t)  -  t  Clai  G(t-T. (,1))  -  •  <  t  <+  . 

J--  •  j  J 

ou  t.  est  I 'instant  os  I'4v4nement  j  ou  processus  oe 
Poisson. 

Les  lois  sph4r iquement  invariantes,  melanges  oe 
lois  gauss'ennes  centre##,  sont  assez  general##  et 
ont  4te  souvent  utilise##  pour  aes  oescriptions  oe 
bruits  marins. 

La  caracterlsation  oe  tel  las  lois,  est  foumie  cans 
to  §  I  4  I'aioe  o'un  theorem#  oe  representation  oe  Bar 
stein  ([5]  p  292). 

Oans  une  seconoe  par  tie,  la  loi  tempore)  le  ou  pro¬ 
cessus  observe  Z  est  oeterminee  par  la  fonction 
caracterlstique.  A  partir  oe  celle-cl  on  montre  oans 
la  Seme  peril#  qua  les  moments  o'orare  impair  sont 
nuls  et  on  calcule  les  moments  o'orare  aeux  et  quatre. 

A  I'aioe  oe  ces  moments,  il  est  possible  oe  calculer 
la  reponse  Impulslonnel le  G  ou  filtre  llneaire,  causal 
e  phase  mlnimale  . 

La  loi  u  oecrivant  la  variable  spheriquement  inva 
riant#  C,  est  approche#  par  une  comblnalson  convex# 


Dirac  u  •  I 


.  ,  “W  a 
k»1  k 


Le  parametre  \  oe  la  loi  oe  Poisson  ainsi  que  les 
2  e  param4tres  et  sont  loentifitfs  oans  la  oer — 
ni4re  partle.  On  utilise  les  moments  E(2?0)^  ) 
p*  1,2, ...,2e  au  processus  7,  statlonnaire  au  sens 
strict. 

Ce  travail  a  <t<  subventlonne  par  oeux  controls 
N  00014-75  -C-  0491  et  N  00014-01-  K  -0373.  Les 
demonstrations  completes  font  parties  o'un  article 
avec  Gualtierotti  Antonio,  article  4  parartre  en  anglais 
Je  tie  ns  4  remercier  Monsieur  C.  R.  Baker,  profes- 
seur  4  I 'University  oe  Chapel  Hill  en  Caroline  ou  More 
pour  ses  nombreusei  suggestions. 


•  -  Caracterlsation  et  examples  oe  variables  aieatoires 

oe  Ipls  spheriquement  Invariantes. 

On  se  llmite  icl  au  cas  ou  I'espace  vector iel  engenore 

par  un  processus  spheriquement  Invariant  est  oe  oimen 

slon  Infinie  ;  en  ce  cas  (cf  [4]  p21),  les  points  extre- 

maux  oe  I'ensemble  oe  telle#  lots  sont  gaussiens,  oe 

oensite  2 

op  -  * 

j  2  j 

0-1)  (x)  -  - — e  ou  ax  est  la  mesu- 

®<  0  /  2" 

re  oe  Lebesgue  .  De  fat; on  naturelle  on  se  restraint  4 
la  classe  oe  lois  spheriquement  invariantes,  oefinle 
4  I'aioe  o'une  probability  u  sur  au  moyen  oe  la 
fonction  caracterlstique  V  u  €  R 


2  2 
a  u 


0-2)  $(u)  -  I  e  2  u(oa) 

[o,-t 

On  peut  appeler  ces  lois,  lois  spheriquement  invarian¬ 
tes  asymptotiques.  Ainsi  ®(o)  -  u  ([p,*!)  -  1  ;  Si  I 'on 
suppose  que  C  variable  aieatoire  reelle  oe  fonc¬ 
tion  caracterlstique  $  est  oe  carry  integrable  alors 


E  (C)  -  0  et 
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( 1-3)  E(C2)  -  ||  C  i!2  -  f  a2  o  u(a)  <  - 

Prenons  pour  probabi  llty  u  ,  une  comblnalson 
convexe  oe  a  masuras  as  Dirac  an  o^,  s2, ...  c#CR+ 
aoit 

(1-41  u  ■u,St.  +  M2lfi*...+u|SJ-e  oil 

at  u,  ♦  U2  +  • . .  +  Ue  -  >  . 

Alnsl  la  loi  aph<r iquement  invarianta  aaaocita  ast 
qyfinie  par  la  fonctlon  caractyristique  : 


2  2 
u  a  . 


2  2 


(1-5)  tp(u)«  U, 


+  p2e  +. . .  +  a 


.  a 

Far  !•  changament  oe  variable  a-4|ic||  »  00 


suppose 

(,“*)  h£iaUa) 

(l-7)<p(u)  -  t 


a  u  (a) 


Las  lols  sph4r  iquement  invariantes  asymptotiques 
ont  une  structure  parti culitrement  simple  qua  nous 
allons  pr<cisar  maintenant.  Rappel  on  s  la  thtforfeme 
oe  Bernstein  (cf  [5]  p  292) 

II  exists  une  masura  positive  v  sur  R+,  telle  que 
V  v  >  o 


iei  PV 

a  o  v  ( b)  si  at  saulamant  si 

* 

♦  ast  k  valeurs  r<allas,  ♦  ast  3  (]o,  *  [)  at 
Vn(H  (-l)n  ~ >  0 


On  alt  que  I  ast  complement  monotone. 
En  effectuant  la  changement  oe  variables 


2 

b«^,  bijectif  sur  R+  on  obtiant  : 


Th<orfcme  I- 1 

Une  loi  oe  probability  P  ast  sphyriquement  inva¬ 
riants  (asymptotique)  sur  (R,3(*  )  si  at  seulement  si 
1*)  P  ast  symytrique. 

2*)tp  est  Inoyflniment  oyrivable  sur  R  \{o) 

2 

3°)  La  fonctlon  f  oyfinle  par  ♦  (u  )»  »p(  v  )  vyrifie 
n 


(-1) 


a  ♦ 


(v)  a  o  V  n  €  Im  V  v€  To,  •  [ 


a  v 

On  vyrifie  que  la  loi  oe  Stuoant  &  k  oegr y>  oe  li¬ 
berty,  est  sphyr iquement  invariante,  asymptotique. 

II  en  est  oe  mime  pour  la  loi  exponential le,  oe  oensite 

f(x)  “  “  e  *  I*  I  avec  X  €R 

2  2  +  2 

En  ce  cas  tp  (u)  -  — 1  at  ♦  (v)  « 

.2.2  ,2 

X  +u  X  +  v 

Ainsi  les  combinaisons  con  vexes  oe  oenslty*  gaus- 
slennas  cantryas  conaulsent  4  aes  Density*  prysentant 
yvantua I  lament  un  point  angulaux  en  x  »  o. 

La  classe  aes  lois  sphyriquement  invariantes  asymp 
totlques  parait  assez  vaste  pour  oycrire  les  varia¬ 
bles  aiyatoires  symytrlques  sur  R  ,  cedes  provenant 
oe  bruit  cw  fond. 

Etuaions  maintenant  Involution  au  court  ou  temps, 
pr4c)sons  la  mooyiisation  ou  processus. 


II  -  Fonctlon  caractyrlstlaue  o'un  processus  oe  Pois¬ 
son  flltry 

Soil  (Tj  >  J  €  Z),  les  instants  aiyatoires  o'un  pro¬ 
cessus  oe  poisson.  Ainsi  pour  tout  interval  la  [a,  b]oe 
R,  la  nombre  aiyatoire  N  a'yvynements  sur  [a,  bi 
vyrifie  pour  tout  n  €  N 

(ll-l)  P  (N  -  n)ta  "Mb'a)  ^  lb~*l  3 -  ou  X 

n  i 

ast  la  paramltre  ou  processus  oe  Poisson. 

L<4vynement  an  T  proouit,  sur  un  capteur,  h  I 'ins¬ 


tant  t,  un  effet  proportlonnel  ft  la  fonctlon  continue 
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G(l,  t).  En  supposant  la  stationnaritt  G(t,  T)  -  G(t  -  T) 
at  la  causalitt  G(t  -  t)  •  o  pour  t  <  T  , 

On  auppoae  oa  plus  qua 

■  .m 

(11-2)  j  lG( a)|2os-/  |  G(t-T)  |2  ot  -  i. 

O  ‘  T 

L'effet  oe  I  '<v<nem«nl  en  t  ,  sur  la  captaur,  aat 
o'tnergle  finie  at  on  a  normalise  G. 

L’tnergie  capita  aat  altatoire  at  la  mookie  consl- 
okre  un  coefficient  C  altatoire  oont  la  lol  aat  rea- 
treinta  4  ttre  ephtriquement  invariants,  ainal  l'effet 
aur  la  captaur  vaut  C  G(t-  t)  avec  G  fonctlon  cer- 
tainc  at  C  coefficient  altatoire. 

Finalamant  on  auppoae  la  llntaritt  oaa  effata  ora 
olvara  tvtnementa  aur  la  captaur  at  I 'observation  Z 
s'tcrit  V  i  {H  • 

Z(t,ui)-  r  C  (x)  G(  t-r  (iu)  ). 

j€*  J  J 

Lea  amplliuoea  altatoires  (C.,  j  €  Z  )  aont  inotpen- 
oantea  oaux  k  oaux  at  lnatpenoantas  ou  procaasua  oa 
Polaaon  (  ;  j  €  ft).  On  montre  assez  facllement  qua 

al 

E(  |  C|)  <  •  at  al  |  |  G(a)|  oa  <  •  alora  Z(l)  aat 

o 

Inttgrable  at 

r  * 

E(  |Z(t)|  )  s  X  E  ( |  C  j )  j  )  G(a)|  oa 
o 

Ca  aont  plut6t  oaa  proprittts  tnergttlquae  qul 
aont  Inttressantea  at  an  raiaon  oa  la  moot llext Ion 
oa  i'obaervatlon  Z  ,  H  aat  poaaibla  oa  otfinlr  entikre- 
mant  la  to!  temporelle  ou  procaaaua  (Z(t)  ;  t  €  R). 
Enon^ons  la  rtsuitat  fonoamantal. 

Thtortme  II-  I 

Si  2  aat  la  procaaaua  oa  Polaaon  filtrt  V  t  €R 
(II  -3)  Z  (t)  -  I  C  G(i-t  )  ai 

r  J  J 

(II-*)  |o(a)  |  «<•  at  al 

(11-3)  I  a*  o  m  (a)  <  *  alora 


Vm€K  «m€" 

V  u-(U)1  u2....um)€  Rm 

(11-6)  $  (u)  «  E(exp  i  <  u,  Z  >  )  • 

axp  -  X  [  oa  o  u(a) 

JR«R 

+ 

C  I  -  axp  t*  2  <  t  uk  G(,k+,>  )2*  1 

ou  <u,  7>  «  u  Z(t  )  +  u„  Z  (t  )+...  +  u  Z  (t  ) 
112  2  fn  m 

L 'ex ten* ion  oa  caa  rtsultats  au  caa  vactoriel  ne 
poaa  aucun  problkme  thtorique  :  I'obaervatlon  peut 

•  tre  vectorial  le  at  provenir  oa  q  capteura.  En  ce 
caa,  la  fonct  Ion  G  eat  vectorial  la  *  q-compoaantaa. 

Ill-  Calcul  oaa  momenta  at  Z 

La  oonnte  oe  la  fonctlon  caracttristique  <p 

permet  la  calcul  oaa  olvara  moments  sous  rtserve 

o'exlatanca  oa  ceux-ci.  D'ailleura  tp  ttant  une  fonc- 
2 

tion  Oe  u  ,  las  moments  o'orare  impair,  aont  nuis. 
L’expression  oa  la  seconoe  fonction  caracttristique 

♦  aat  particular  ament  simple  : 

(III— I)  f(u)  -  log<p(u)  ■  -  X  f  oa  <%i(a) 

JR»R, 

[l-axp-'Y  <u,G>2(s)] 

ou  I 'on  a  post 

(III— 2)<ru, G  >  (a)  •  u,  G(t+s)  ♦. ..+  u  G(t  +s). 

It  mm 

Las  otve loppemants  an  stria  oa  la  fonction  «p  an 

fonctlon  oaa  moments  pula  oa  I 'exponential  le  dans  la 

seconoe  fonction  caracltristlque  conouisent,  en  tga- 
2 

tant  lea  puissances  oe  u,  at  oe  u.u„,  4  : 

I  I  2 

Proposition  lit—  I 

(111-3)  E(Z(t)2)  -  X  [  a2  ou(a)  [  G(s)2  os  V  t  €  R 

Jr 

(III-*)  E(Z(t))Z(ta))- xj„  a2ou(a)JR  G(t  |+s)G(t^s)os 
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Oi  rappelle  que  I 'on  a  pos t  G(»)  as  -  1  et  que 

JR 

G(s)  -  G(s)Ks)  .  On  voil  que  pour  tout  t  »  t  la 

to.-[  '  2 

formule  (4)  oonne  (3). 

En  effectuant  les  o4veloppements  Ilmlt4s  jusqu'4 

4  2  2 

I'orore  4,  nous  4cr Ivons  l'6galit£  oe  u  *l»u.u. 

112 

Proposition  I II  —  2 

E(7(t)*|  •  )  l  f  a*<si(a)  f  as  G4(s)  + 

JR  +  JR 

3  X  2  [j  a2  Oi(a)  f  as  G2(sl  ]2 
JR  +  JR 

E(Z<t])2  2<t2)2  )  - 


3X  jR  a4ou(a)jRas  G2(t]+  s)  G2(t2+s)  + 


2 

+  X2^J  a2  csa(a)  j  as  GZ(s)  j  + 


+  2X”[f  a“  Oi(a)  !  asGlt+sl  Gd.tsl 
L*R+  JR  '  2 


IV  .  loentlfication  oe  la  fonction  ae  Green  G 

La  covariance  T  ae  Z(t),  est  stntionnaire  et  vaut 

V  t  6  R  : 

(IV- I)  r(T)  -  E(Z(t)2(t-r)  )  » 

X  JR  a2  o  u  (a)  J"  G(s)G(s+t )  as 

f  2  +  ^ 

Ainsi  \ ,  bom  (a)  est  un  facteur  oe  normalisation  ou 
”R+ 

X  est  le  paramitre  oe  la  toi  oe  Poisson  et  ou 

a2o  u(a)  est  le  moment  o'orare  2  oe  I'amplituoe 

JR+ 

allatoire  C,  sph4riquement  invariante.  Comme  I 'on 

a  suppose  j  |  G(s)  |  2  os  -  I  ,  on  a 
r  *+2 

T(o)  ”  X  J  a  a  u  (a). 

La  fonction  G  4  identifier  v4rifie  oonc  I'4gallt6 


(IV—  2)  *  -  [  G(s)  G(si-T)  as. 

r<o>  jr 

La  transformie  oe  Fourier  va  nous  Conner 
un  moyen  oe  caiculer  G  sous  reserve  a'hypothfeses 
compl6mentaires. 

Proposition  'V-  I 

+• 

Soit  g(f)  «  3  (G)  (f)  ■  |  e,afG(s)as.  Alors 
(IV-3) 

r(T)-rto)  lm  e-l*'|g(f)|2  *£• 

Ainsi  oans  le  langage  oes  automatlciens,  le  filtre 
causal  oe  r4ponse  impulsionnel le  G  ,  est  oe  galn|g| 
connu  :  il  suffit  a'effectuer  la  transform4e  ae  Fourier 
inverse  sur  V  .  La  oom6e  ou  gain  |g|  ne  determine 
cepenaant  pas  la  rlponse  impulsionnelle  G. 

Definition  IV-2 

Un  filtre  causal  G  est  4  phase  minimale  si  1/g  est 
4galement  la  r6ponse  fr<quentiel le  a'un  filtre  causal. 

En  ce  cas  le  gain  |g|  determine  le  filtre  causal  4 
phase  mlnimale  (cf  [3]  p  36) 

Pour  g  rlponse  frequent lelle  a'un  filtre  causal  4 
phase  minlmale  une  extension  analytlque  oans  R+  est 
possible  et  on  o4finit  V  b  >  o 

_+  m 

GbU)  “  e  Ifl  g(f+  lb)  o  f  . 

b 

Or  cette  fonction  G  s'exprime  uniquement  en  fonc¬ 
tion  ou  gain  j  g  jpar  ('expression 
(IV- 4)  Gb(t)  - 

_L  f  “ift  I  f  Y  (f+ib)+l  loq(q(y)  i 

2tt  J  016  **P  ini  ,  ,  J  2  '  °  > 

~m  R  Y-(f+'b)  Y  +1 

Pour  b  assez  petit,  nous  avons  ainsi  une  approxi- 
b 

mation  G  oe  G  ,  rlponse  impulsionnelle  o'un  filtre 

Iln4aire,  causal  4  phase  minlmale,  en  effet  Gb  conver- 
2 

9*  o«n*  L.  (R,  t  at)  v*r*  G  torsque  b  tcna  vers  O. 


IDENTIFICATION  a'UNE  CLASSE  DE  PROCESSUS  DE  POISSON  FILTRES 


V-  loentification  oe  X  et  oe  u_. 

La  fonction  G  est  suppos^e  coonue,  elle  a  iti  ioen- 

tlfiye  aan s  le  paragraphe  ant^rieur  et  ainsi,  pour  tout 

t  2p 

p  6  K  ,  on  con-taft  la  valeur  num^rique  oe  G(a)  as 

JR 

De  plus  le  oyveloppement  oe 
♦  (u)  •  log  ®(u)  - 


I  u2*5  UJ2_  X  [  a2**  o  u  (a)  j  G(s)2p  os 
p« t  2^p!  B+  R 


cooouit  4  poser 


(V-  I )  X  a  ^ 0  u  (a)  »  cr  p*l,2 .  2e 

JR+  p 

L'estimation  oes  moments  E(7(t)2*>)  pour 

p  -  1,2,...,  2e  fournit  le  a^veloppement  oe  cp  et 

4e 

ensuite  oe  t  jusqu'au  terme  en  u 


(V-4)  P^l-a)  -  P2k(a) 


(V/-5)P2k+,(-a) 


P_.  .<«> 

2k+  I 


Disposant  oes  seules  constantes  afJt  or  ,  .  .  .  ,  a  , 
on  cherche  une  approximation  ia  oe  la  forme 
(V-6)  u  «  u,  +  Uj  S  0  2  +. . .  +  ue  ^0^ 
comblnalson  eonvexe  4  a  termes  ,  oe  Dirac,  situ£s 

Bn  °l*  °2 .  °e  «  V 

On  suppose  lei,  qua  I'amplltuoe  al^atoire  C  a  pour 

loljle  melange  oe  e  lois  gaussiennes  centres,  est 

oonc  sphifrlquement  Invariante  oe  fonction  caract^ris- 

2  2  2  2 
tique  :  u  a  .  u  o 

_ ]_  e 

_  ,  I  u  C,  2  2 

E  (e  )  -U,  +  ...  +  U. 

Les2e+I  inconnues  Wj,...|i  ,  e*  *  v4- 


Dans  ce  paragraphe  ,  il  s'agit  oonc  a'approcher  la  riflent  las  2e<-l  Equations  o^ouites  oes  £galit£s  (V-l) 
probability  a  ytaiye  sur  R+  ainsi  que  le  paramfetre  X  at  (V-8)  : 
oe  la  loi  oe  Poisson  4  partir  oes  <galltys  (V- 1)  l^l+  ^2  +  *  "  + 


Si  le  paramitre  X  oe  la  loi  oe  Poisson  <tait  connu, 
le  probifeme  se  ram4neralt  4  une  situation  classlque  : 
Introouisons  la  probability  u  symytrlque  sur  (R  ,  X  ) 
et  oont  la  restriction  u/  R+  4  (R+1  l5l+)  soit  pr4cisy- 
mcnt  u  .  De  cette  fa^on,  les  moments  oe  £  sont 
connus  : 

(V-2)  J  a2**  unit)*  j  a^o  u  (a)  «  /  X 

(V-3)  [  a2***'  ou  (a)  -  0. 

JR 

L'ioentification  a'une  probability  u  4  partlr  oe  ses 
moments  est  un  probl4me  rysolu  (cf  [I]  )  et  nous 
particularisons  les  rysultats  g4n4raux  au  cas  a'une 
probability  svmytrigue  u  sur  (R  ,  X  ). 

L'orthonormalisation  oe  Gram-Schmia*  °*  la  suite 
2 

(I,  a,  e  , . . .  )  concsjit  aux  polynflmes  orthogonaux 

<p0(a>  ,  Pt  (a),  p2  (a),...)  vyrifiant 


X  [U.°t  +^2ff2  + 


...  Ue02J-  ■>, 


Ur  »  2. 2e  ^  _  2.  2e^  2.  M 

|M*t9l  +  *2*2  +-+Vt1  =*« 

L'ytuae  oe  ce  syst4me  non  linyaire  conouit  aux 
rysultats  suivants  : 

X  est  solution  oe  liquation  linyaire 


(  V-8)X 


Q_  .  .  . 

or 

2 

e+  1 

O 

m  • 

e-f  I 

2e 

O 

a,  ... 

Of 

e 

Q, 

Qf  ... 

Qf 

1 

2 

e+  1 

a 

Or  .  .  .  . 

O  i 

e 

•♦1 

2e  , 

Cette  formule  est  originals. 

Maintenant  que  X  est  connu,  on  applique  les  rysul¬ 
tats  class!  ques  ou  probl4me  oes  moments. 

Le  pol  ynome  P  ,  oe  oegry  e  en  u2  ,  a  pour  zyros 

cans  R  pr£ci  foment  C  .  ,  0_ . a 

*  I  2  a 


IDENTIFICATION  D'UNE  CLASSE  DE  PROCESSUS  DE  POISSON  FILTRES 


Avec  X  et  connus,  le  systime  (V-7)  est  Iin6aire 

en  u.  avec  reoonoance  oes  Equations.  On  le  r6souo 
k 

ou  bien  I'on  utilise  les  ygalit6s  (cf  [l]  p  22  ). 


j_  _ _ !_ 

2  <rl  k  *2  e-  I 


1 1  2y  •  •  •  ,® 


En  conclusion  ,  la  moctyiisation  Introouite  aans  ce 
travail,  sous  reserve  o"hypothises  assez  natureltes, 
permet  oe  o6finir  entifcrement  la  loi  temporal  le  a'un 
processus  observation  Z  (vector idle  ou  non).  Les 
transformations  non  Iin6aires  sur  7  peuvent  oonc 
#tre  <tuoi<es.  Nous  sortons  ou  caore dee  processus 
ou  seccono  orare  et  oes  seules  oenii  t6s  spectrales 
oe  puissance.  L'hypothtse  gaussienne  peut  8tre 
abanoonn6e  . 

Le  moafcle  anttfrieur  permet  o'ajuster  toute  cova¬ 
riance 
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T(t)  »  E(Z(t)  Z  (t-  T)  )  t  €  B  et  tout  moment 

2p 

E(7  (t)  )  p  -  1,2, ...  ,  2e  ,  les  moments  Impairs 

2k*  I 

E(7(t)  );  k  €  h  <tant  supposes  nuts. 
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